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EMISSION, ABSORPTION AND POLARIZATION OF GYROSYNCHROTRON 
RADIATION OF MILDLY RELATIVISTIC PARTICLES 


ABSTRACT 

In this paper, we present approximate analytic expressions for the 
emissrvxty and absorption coefficient of synchrotron radiation of mildly 
relativistic particles with an arbitrary energy spectrum and pitch angle 
distribution. From these, we derive an expression for the degree of 
polarisation. To accomplish this, vc use methods of integration developed in 
a previous paper. 

Ve tnen compare tne analytic results with numerical results for both 
thermal and non-thermal (power law) distributions of particles. 



I. INTRODUCTIOH 


The forniilaa for evaluation of emiativity and absorption coefficient of 

synchrotron radiation in the ultra-relativistic (synchrotron) and non- 

relativistic (cyclotron) limits have been known for decades. * • In the 

intermediate energy range, however, no simple formula exists for an arbitrary 

distribution of particles. Calculations have been made of the absorption 

3 4 

coefficient of a thermal plasma * , but these results are limited and not very 
convenient for quantitative calculations. As t result, the usual practice has 
been to use lengthy numerical calculations.^*^. 

In general, the low harmonics of cyclotron radiation will be self- 
absorbed, absorbed by the surrounding plasma, or suppressed by the Rasin- 
Tsytovich effect. Consequently, we are interested only in radiation at the 
higher harmonics, where the optical depth is less than or equal to 1. 

In a recent paper^ we presented simple approximate methods for evaluation 
of the frequency spectrum and dependence on observation angle of the 
synchrotron radiation at high harmonics from an (essentially) arbitrary 
distribution of particles in a given magnetic field. In this paper, we shall 
use the same methods to calculate the emissivity and absorption coefficient of 
the extraordinary and ordinary modes of synchrotron radiation separately for 
an arbitrary particle distribution. Also, we will derive the degree of 
circular polarization from these expressions. (We will refer to the paper in 
which these methods were developed as Paper I.) 

In tne next section we describe the general results and in Section III we 
use them to find emissivity, absorption, and polarization of radiation from 
particles with Maxwellian and power law energy distributions. Section IV 
contains a final summary. 
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II . GENERAL RESULTS 


Consider particles with charge e, mass and distribution f(M*Y) where 
fdpdy is tne number density of particles in the energy intervals in the energy 

n 

interval (in units of m^c^) from y bo Y ^ with pitch angle cosine 
between y and y + dy. The emissivity and absorption coefficients for the 
ordinary (+) and extraordinary (-) modes are^» 



where 


n. 


1 “ ( X j'(x ) 

— 5 — / I a+(cos6-3li)J (x ) - (1 - pycosB) ---- — 

„ 2^, 1 i mm 

*+ m=l I ® 


6(y). (2) 


-3y* 9 /f(P.Y)' 


3cos0-y 


1 9f(y,Y) 


0) 


f(y,Y) 9 y V 3y^ / Y3^(i-y*)^^^ f(y.Y) 9y 


(3) 


mVb 

y B _^(1 - 3ycos0), 3^ = (vVv^)3sin6(l V^, “ eB/2mgC . (4) 

The quantities represent th,. ratio of the semi-major to the semi- 
minor axis of the polarization ellipse for the ordinary and extraordinary 
modes and, in general, are complicated functions of angle 9, v, and the 

plasma frequency V^. However, when the frequency V >>'Vp (consistent with 
<1), these simplify to 

a+ “ x/fl +(1 + X « -2vcos0/(Vj^sin20) , (5) 
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Since we are only dealing vitn the higher harmonica, b is large and we 
can replace the sum in equation (2) with an integral. This integral can then 
be done using the 6- function. Also, since m is large, we can approxiaate 

Q 

the Bessel's function as . 






Z - 


ze 


(l-z2) 


1/2 


(2inn) 


1/2 




( 6 ) 


which is valid for m(l - » 1. 

With this approximation the integration over the pitch angle is carried 
out using the method of steepest descent. This gives excellent results as 
long as the pitch angle distribution is not extremely anisotropic; 
dln£(p,Y)/dy >> (See paper I, eqs. 7 and 19). This also enables us to 
drop the second term in the expression for O) in equation (3) and amounts to a 
substitution setting y 3 cosO . 

We also use the method of steepest descent for the integral over the 
energy. Evidently, this is not as well justified as using steepest descent 
for pitch angle, but it does give an excellent approximation for high 
harmonics and particle distribution which fall rapidly with increasing energy, 
(e.g. a power law foc^y-1) ©r a thermal distribution f«e ' 

Doing these integrations, we find for j^ and ; 


1 

2 


V 

l.~p- (V/Vj,)'/' Y,(l^„co. 
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r . .L r 1 . 

a* + 1 L 


cot 0 2a 

+ 1) (l-$* cose.) 

yt y 

0 'o 


1/2 
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( 8 ) 


X‘2 . Y* 

® dy" 


dln C 
Y„ dlny 


W(t„.6) 


(9) 


® (1 + t^^)» 

Vb 


l/(l+t 

t^e o 

o 


( 10 ) 


Note that all these expressions are evaluated at critical energy (and 
the corresponding 0^ and t^ ** Yq8q sin 9). These are two^ such critical 
energies; one for emission and one for absorption. These are obtained from 
the transcendental equation; 

^These critical values are the same for both modes since in our approximation 

the parts of the expression which are different for each mode are slowly 

varying functions of y and y and do not affect the values of the critical 

parameters u 3 cos6 and y 
0 0 o 


t^(l + t^)"^^^+(l - e^co8^0/t^)lnZ^jj “ 
”(Vb/2v sin^e )(dlnC/dy) = e/sin^e . 


( 11 ) 


Here and in equation (7) the function and are 

related to the particle distribution f as 

/ f($^ cos9,y)\ 

C • » f (0 cos0,y)/y , C “-0yd/dy( I (12) 

JO K \ 0y* / 
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In tbe expression (eq. 9) for IC , the function W is a complicated 
function of t^ and 6 . However, as shown in paper I, W can be 
approximated by 

W - 3/2 + 1/(y* - 1). (13) 

0 

This is an excellent approximation in the entire region where equation" 
(61 to (ill are valid. 

Polarization . For a single particle, the synchrotron radiation can be 
represented by the Stokes Parameters^; 


TTj - + n_ , 

iTq -[(1 - a^^)/(l + a+^)](ri^ - n_) + cos6[a+/(l+a^^) ], 

• 2Vn_^n_ sin6 , 

- l2a^/a+a^^}l(n^ - n_) - 2^Jri^_ cos6[(l-a+^)/(l+a+^)] , 


(14) 


where and a+ are defined as before, and 6 is tbe phase difference 

between tne modes. 

In tne limit of large Faraday rotation (negligible absorption over a path 
length in which the plane of polarization rotates by 2ir), we have the time 
average values <cos5> “ <sin5> “ 0 . 

Thus n ** 0 
u 


_ 1 - a; - 1 

“ Ci( n - n ) , c. - 2 *“ 2 1/2 


-X 


n„- 2c2(n. - n_), c„ = 2 “ 

^ ^ ^ 2 1 + a+ 2(1 + 


(15) 
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From this ve obtain the degree of polarization 

p . . K . 

Hj + nj 


( 16 ) 


For a distribution of particle ve integrate over the distribution as 
before and find P ■ Jj^ - j_| /(j^. + j«) which with the help of equation (7) 
reduces to 

[1 - 2xcos6(l - 3^9co8^6)^^^/Y^sin^6] 

P 0 ^ 

(1 + ♦ 2 cot^e/Y' ) 

Equations (7) to (13) along with (17) give our results in their general 
forms and are valid for all particle distributions which are not extremely 
anisotropic. (The results for the extremely anisotropic situation are more 
complicated and were described in Paper I. The modification obtaining the 
emission ana absorption coefficients of the two modes separately and the 
polarization is similar to the modification described above. We will not 
present these results here because of their complexity and their limited 
usefulness.) 

Given a distribution function subject to these limitationsp the first 
step is evaluation of critical energies Yq from equation (11). Then 
equations (b) to (10) and (12) and (13) evaluated at the appropriate Yq's 
give the desired results. This is still a complicated procedurcp especially 
the solution of equation (11) for 

In toe next section ve shall show how this is simplified considerably for 
the two most commonly used particle distributions. Before doing so we 


6 



consider the nsymptotic limits of these equations. 


Asymptotic limits . Let us consider first the case vhen angle 8 is 
not too small (i.e. radiation avay from the direction of the field). Then in 
tne two extreme limiting cases* equation (11) simplifies 


i) e « 1, » 1» B*Yq ■ 2/3e sin0 , W ■ 3/2; 


(18) 


ii) £ » 1, 0Q « 1. Yq If PJyJ - 4/e . W - c /4. 


The first case is realized at high frequencies and particle distributions 
which are not extremely non-*relativistic which is the case of interest here. 
The seond case is valid for non-*relativistic particles and at low frequencies 
and has limited usefulness. 


111. EMISSIVITY AND ABSORPTION COEFFICIENT OF TWO COMMONLY USED 

PARTICLE DISTRIBUTIONS 

The two particle distributions we use as examples are i) the distribution 
from a thermal gas* i.e. a Maxwellian distribution in energy and isotropic 
pitch angle distribution; and ii) the distribution with a power law spectrum 
at high energies and with a slowly varying pitch angle distribution. 

A. Thermal Spectrum 

3 

This is the type of particle distribution considered by Trubnikov . 

A 

The distribution of particles at temperature kT (in units of m^c ) is 
f(p,Y) " Y(Y^ - 1 )^^^ , (19) 

where for kT 1 
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> 


( 20 ) 


C P 0x/2)[2TT(kT)T*^ (1 - 15 kT/8 + ...) 


and n la tne number density of particles* From these and equation (12)y it 
is clear that kTC - C- - £/Y and that dlnC /dy - dlnCi/dy . Thus, it is 

fC J Is J 

obvious from equation (11) that the critical Yq same for both and 

and. in fact, according to eq. (7), “ v*kTK^. We also find 


dln(f/y) 

dY 



kTy 

Y-1 


d^ln(f/Y) 


-(y* ♦ 1) 
(y* - 1)^ * 


( 21 ) 


Using these equations we can calculate Y^ X from equations (11) and 
(9) respectively. As shown in Paper I. these expressions can be considerably 
simplified. We find that the following expressions 


r(2vkT/v^)U + A. 5 vkTsin^/vj,)”^^^ kT < 1 




(4'OkT/3Vj,sin0)2^^ 


kT 5S 1 


( 22 ) 


- (2kT/Y^)( Yq - D/OYq “ 1) 


kT 1 


have tne correct asymptotic limits in agreement with equations (18) and agree 
with the exact results from equations (9) and (11) to within 30Z for most 
relevant ranges of angles, frequencies end temperatures and better than lOZ in 
the majority of the interesting cases. 

The above equations and equations (7), (8), and (10) give a complete 
description of the emissivity and absorption coefficient from a Maxwellian gas 
at all temperatures and frequencies. They are valid for kT < 1 because at 
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temperaturea kX > 1 the use of the method of steepest descent for integration 
over the energy becomes invalid. However, the existence of the extremely 
relativistic thermal gas is in doubt^^. On Figure 1 we compare the total 
absorption coefficient K -t- K_ obtained from these relationships with 
numerical results from Lamb and Masters^. As evident, our analytic results 
give excellent agreement to the detailed numerical results even at low 
harmonics . 

In the two limiting cases described in the previous section, these 
equations are considerably simpler. The interesting case, c « 1 corresponds 
to VkT/Vjj » 1, Y* " 4vkT/3VjjSin6t so that 


j + 


^kTK^ - (2^^^7Te\)j^/3c)C(vkT/vj,) 
- 1 qi2(sine)^^^ (3v^/4vkl)^^^ 



(23) 


B. Power Law Energy Spectrum 

Power law spectra are commonly used spectra in astrophysical problems and 
in other problems when the tail of the Maxwellian distribution begins to 
deviate from the exponential form. Usually power law spectra are defined with 
a low energy cut-off. To avoid such discontinuities and the divergence of the 
number of particles, \'e assume a spectrum of the form 

f(M,Y) - Cg(M)tl + (Y - 1)/ ec)]'^ (24) 


which converge at low energies; 


C “ n(6 - 1/e^) , 


r+1 

I g(y)dy - 1 . 


(25) 
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From this and equation (12) we have 



c^(y)*i<y ) . 


£/Y nnd 


( 26 ) 


In (24) plays the role of the low energy cutoff (in units of m^c^). For 

energies much greater than the spectrum is a power law with index -6 

bur it tends to a constant value at lower energies. The particles can be 

classified as ultra-relativistic or non-relativistic if £ >> 1 or C << 1. 

c c 

We are interested primarily in cases with sk 1. 

For distributions which are not highly anisotropic (i.e. dlng(U)/dV« 
v/V|j)» we can carry out a calculation similar to that for a thermal gas. From 
(24) and (26) we find 


— d InC j 

6y 

^ 1 a. 


Y^d^inCj 

dlnY 

" 1 + 

Y-l+e^ 

$ 

dY* 

and 




-dlnC^ 

-dlnC 

y^>^(y) 

Y*d*lnCj^ 

dlny 

dlnY 

*j(y) 

dY^ 

where 





d^ (y) 


d$ (y) 

4, (Y) - . 

4> (y) 

- ~2 . 

Z 

dY 

2 

dY 


6y* 

1 + (27) 

(y-l+e )^ 

' c 


Y^d^lnC. \ (y)$. (y)-^*(y) 
1 + X i 2 (28) 

dY* ^2(Y) 


(29) 


It turns out that for semi-relativistic particle energies, for 

>> 1 we can use equation (18), and the values of X and Yq obtained from 
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(11) and (9) can be approximated aa follows 


eJVo -[4v/(3v^sln6)l 


! 1/(1 ♦ 5) for j + 

1/(2 ♦ 6) for Kj 


( 30 ) 


and 


(1 + for 

(2 + 6)”^^^ for «c^ 


(31) 


These expressions are valid to within 30Z for ^ ^ ^ ^ 

sin6 «« and v > Vk* However, they are most accurate at high 

frequencies v >> where » 1 and 3^ « 1. In this limit, 

substitution of (30) and (31) into (7) gives 


/ Nr3ireVsin0 
b 


v*k+(v,6) 


Ac 


Ce^g(B^cose)Y^(0) 


[(3e*v (6-f-l)sin0/4v)^^"^^^^ ^-(6+l)/2 

C D 


(32) 




'c b 


\"c / 


In this limit and their difference is of higher order 

in Y,, ; 

2sin0 

Y^(0) - IT . (331 

These expressions have exactly the same dependence on and 6 as 

2 

the ultra-relativistxc forms . 

For extremely non-relativistic particles, that is, for « 1 the above 
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expreaoions are valid as long aa vA^^5 » 1. For the unlikely case of vA>|^6 

« 1, ve find Chat y ta i and e» vc /vk 5> Substitution of this in (7) 
’ 'o c o 

gives tne emissxvity identical to that expected from a thermal gas if one 

identifies vith the temperature kT . 

When ve substitute eqs. (28) to (31) into (7) and (19)» (using 6 <■ 60^, 

d *4), ve get the results for and P seen in Figures 2, 3, and 4. 

8 

These are sbovn in comparison vith numerical results of Marsh and Dulk . Note 
that in this comparison ve have a distribution vith a sharp cut~off point at 
c - 0.02, so tne approximations for y^ at V « are not as good as ve 

vould like. Still, at high v/\^ we get to vithin 5Z in general. 

The values for and are vithin 60Z at high frequencies and the 

approximation is better at higher The polarisation looks better at 

lover frequencies (although it still is better than 50Z even vhen P is small 
and errors can be magnified). The results for the extraordinary mode are not 
as good as those for the ordinary mode. (The percentage of error never 
reaches 30Z for either or above 10^’ V Also our results are 

(nearly) systematically higher than the numerical results. 


IV. SUMMARY 

Using a simple metnod of integration developed previously, ve have 
derived expressions for the emissivity, absorption coefficient and 
polarization of synchrotron radiation for an arlitrary distribution of 
particles. 

Equations (7) to (13) and (17) give our results in their most general 
form. And we find that equations (9) and (11) can be simplified considerably 
as in equations (22), (30), and (31) for the thermal distribution and the 
povcr law. 
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Our results do sgree vith previous snslytic results* end they give good 
epproximations to detailed nuaerical results. Although our results vere 
derived for high hareonics, they give good agreeaent down to lover haraonics 

for V lOu^ for for total j and tc * however* the agreeaeat 
is good down to V 

These results are Halted to pitch angle distributions which are not 
extremely anisotropic and energy spectra that decrease rapidly with 
increasing energy. They also are only applicable for ealssivities and 
absorption coefficients away from magnetic field times. Examples of other 
cases can be found in Paper I. 

A more detailed numerical comparison aill be found in a subsequent paper 
to be published in The Astrophvsical Journal. 
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FIGURE CAPTIONS 


Figure 1. The total ayuchrotron abaorption for a theraal aource at 6 •t/1 
for kT “ 0.04 (20 keV elccCrona). Pointa are fron analytic expreaaiona; the 
aolid Ixnea are nuiaerxcal reaulte of Lanb and Kaatera^. 

Figure 2. Synchrotron eniaaivity of each node divided by nagnetic field and 
total particle number. Log (J^/EN) va Lojt (v/v^)» at 0 •• 60*^» 6*4. 
The ordinary mode haa been ahiftcd down by a factor of 10 for clarity. The 
aolid linea are numerical reaulta of Dulk and Marah. The O'a are our analytic 
reaulta. (The aame will be true in Figures .1 and 4.) 

Figure 3. Self abaorption of each node. Log (Bic. /N) va Log (v/vj. The 

” o 

low pointa at v/v^ ” 10 are due to the fact that we uaed a diflcrent 
approximation for on the first two pointa; one which ia better at low 
v/v,,. 

Figure 4. Degree of Circular Polar iaat ion va l^g (v/v^). 
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